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4.  INTRODUCTION 


The  research  project,  “Time-Resolved  and  Spectroscopic  Three-Dimensional  Optical  Breast 
Tomography,”  is  devoted  to  developing  a  safe  and  affordable  breast  cancer  detection  modality 
that  makes  use  of  noninvasive  near-infrared  (NIR)  light  for  imaging,  locating  and  diagnosing 
tumors  in  human  breast  with  high  resolution  and  specificity.  The  research  focuses  on  the 
development  and  application  of:  (a)  experimental  approaches  for  probing  the  target  using  NIR 
light,  and  (b)  image  reconstruction  and  target  localization  algorithms  for  obtaining  target  location 
and  generating  three-dimensional  (3-D)  tomographic  images. 

Significant  progress  has  been  made  in  both  the  areas  during  the  reporting  period  (May  1,  2007  - 
April  30,  2008)  covered  by  this  report. 

5.  BODY 

The  tasks  performed  and  the  progresses  made  during  the  current  reporting  period  are  as  follows: 

•  Application  and  further  refinement  of  optical  tomographic  imaging  using  independent 
component  analysis  (OPTICA)  developed  during  the  earlier  reporting  periods  for  locating 
and  cross-section  imaging  of  a  tumor  in  a  model  breast  constructed  using  ex  vivo  breast 
tissue  specimens; 

•  Development  of  a  new  approach  for  detection  of  a  target  in  a  highly  scattering  turbid 
medium  using  near-infrared  center  of  intensity  time  gated  imaging. 

We  provide  a  brief  outline  of  our  accomplishments  in  these  areas,  and  refer  to  appended 
publications  for  detailed  description  where  applicable. 

5.1.  Detection,  localization,  and  cross  section  imaging  of  model  cancerous  breast  using 
OPTICA 

We  have  pursued  further  improvement  and  testing  of  OPTICA  (optical  tomographic 
imaging  using  independent  component  analysis)  that  was  developed  ’  and  reported  during  the 
earlier  reporting  periods.3,4  Particular  emphasis  was  placed  on  testing  the  efficacy  of  the 
experimental  arrangement,  theoretical  model,  and  numerical  algorithm  to  detect,  locate,  and 
image  a  tumor  inside  a  cancerous  breast  (TO  1,  Task  2;  TO  2,  Task  4;  and  TO  4,  Task  9).  These 
are  detailed  in  Appendix  1:  “Optical  diffuse  imaging  of  an  ex  vivo  model  cancerous  human 
breast  using  independent  component  analysis,”  IEEE  J.  Select.  Topics  Quantum  Electron.  14,  43 
(2008). 

At  our  present  state  of  development  OPTICA  can  detect  an  optical  inhomogeneity  (such  as,  a 
tumor)  inside  a  highly  scattering  turbid  medium  (such  as,  a  human  breast)  based  on  a  difference 
in  optical  absorption,  fluorescence,  or  scattering  characteristics.  ’  ’  The  multi-source  probing  and 
multi-detector  monitoring  scheme  used  in  OPTICA  provides  multiple  angular  views  of  the  target 
and  enable  determination  of  the  position  of  the  target  with  high  accuracy  (typically  within  1 
mm).  The  fonnalism  has  been  further  extended  to  provide  cross  sectional  image  of  a  target 
within  a  turbid  medium. 5,6 

We  will  provide  a  brief  overview  of  a  measurement  on  a  realistic  cancerous  breast  model  and 
leave  the  details  to  a  recent  appended  reprint  {Appendix  1).  The  model  breast  was  a  70  mm  X  55 
mm  X  33  mm  slab  composed  of  two  pieces  of  ex  vivo  human  breast  tissues  provided  to  us  by 
National  Disease  Research  Interchange  under  an  Internal  Review  Board  approval  at  the  City 
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College  of  New  York.  The  larger  piece  was  normal  tissue  that  included  mainly  adipose  tissue 
and  streaks  of  flbro-glandular  tissues.  The  existence  of  the  fibro-glandular  tissues  was  not  known 
prior  to  making  the  measurements.  The  second  piece  was  mainly  a  tumor  (infiltrating  ductal 
carcinoma)  with  a  small  amount  of  normal  tissues  in  the  margins  with  an  overall  approximate 
dimension  of  8  mm  X  5  mm  X  3  mm.  An  incision  was  made  in  the  mid-plane  (along  the  z-axis, 
which  was  the  shorter  dimension  of  the  tissue)  of  the  nonnal  piece  and  some  amount  of  normal 
tissue  was  removed  from  the  central  region  making  a  small  pouch.  The  tumor  piece  was  then 
inserted  into  the  pouch  and  the  incision  was  closed  by  moderate  compression  of  the  composite 
consisting  of  the  normal  tissue  and  the  tumor  along  x-y-z  directions.  The  breast  tissue  slab  was 
contained  inside  a  transparent  plastic  box.  One  of  the  sides  of  the  box  could  be  moved  to 
uniformly  compress  the  tissue  along  the  z-axis  and  hold  it  in  position.  The  resulting  specimen,  a 
70  nun  X  55  mm  X  33  nun  slab,  was  treated  as  one  entity  in  the  subsequent  imaging  experiment. 
The  position  of  the  tumor  within  the  slab  was  known  since  it  was  placed  in  position  as  discussed 
above.  One  of  the  tests  of  the  efficacy  of  this  imaging  approach  was  to  see  how  well  the  known 
position  is  assessed. 

The  experimental  arrangement  (shown  schematically  in  Fig.  1(a)  of  Appendix  1 )  used  a  200- 
pm  optical  fiber  to  deliver  a  784-nm,  300  mW  continuous- wave  beam  from  a  diode  laser  for 
sample  illumination.  The  beam  was  collimated  to  a  1-inm  spot  onto  the  entrance  face  (the 
‘source  plane’)  of  the  slab  sample.  Multiple  source  illumination  was  realized  in  practice  by  step 
scanning  the  slab  sample  across  the  laser  beam  in  a  22X  16  x-y  array  of  grid  points  with  a  step 
size  of  2.0  mm  using  a  computer  controlled  translation  stage.  The  signal  from  the  opposite  face 
of  the  sample  (the  ‘detection  plane’)  was  collected  by  a  camera  lens  and  projected  onto  the 
sensing  element  of  a  cooled  16-bit,  1024  X  1024-pixel  charged  couple  device  (CCD)  camera. 
Although  the  scanned  area  is  42  mm  X  30  mm  on  the  source  plane,  the  imaged  area  of  the 
detection  plane  was  much  larger,  covering  the  entire  70  mm  X  55  mm  area  of  the  model  breast. 
Each  illuminated  pixel  of  the  CCD  camera  could  be  regarded  as  a  detector. 

For  illumination  of  each  scanned  point,  the  CCD  camera  recorded  an  image.  A  typical  raw 
image  is  shown  in  Fig.  1(c)  of  Appendix  1.  Each  raw  image  was  then  cropped  to  select  out  the 
information-rich  region,  and  binned  to  enhance  the  signal- to-noise  ratio.  All  the  binned  images 
corresponding  to  illumination  of  the  grid  points  in  sequence  were  then  stacked,  and  used  as  input 
for  independent  component  analysis.  The  details  of  the  analysis  method,  theoretical  formalism, 
target  localization  algorithm,  and  experimental  arrangement  have  been  published,6  and  are 
presented  in  Appendix  1.  After  optical  measurements,  the  sample  was  transferred  to  our 
collaborators  at  the  New  York  Eye  and  Ear  Infirmary  for  pathological  study  and  correlation 
{TO 2,  Task  5). 

The  key  results  of  the  experiment  are  as  follows. 

(a)  OPTICA  identified  three  different  structures  (Fig.  2,  Appendix  1 )  as  three  independent 
components  based  on  the  differences  in  optical  properties.  These  include  the  tumor  whose 
presence  and  position  was  known  from  the  sample  preparation  process.  We  ascribed  the  other 
two  structures  to  fibro-glandular  tissues,  as  the  normal  component  of  the  breast  tissue  specimen 
was  mainly  adipose  tissue.  Comparison  with  the  pathology  results  further  confirmed  the  identity 
of  the  tumor  and  the  fibro-glandular  tissues. 
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(b)  The  location  of  the  tumor  was  determined  to  within  ~1  mm  in  all  three  dimensions.  The 
locations  of  the  fibro-glandular  tissues  were  also  estimated.  The  locations  of  the  components  are 
given  in  Table  I  of  Appendix  1. 

(c)  The  FWHM  of  the  tumor  is  estimated  to  be  ~  10.3  mm  and  7.4  mm  along  the  x  and  y 
directions,  respectively  (details  in  Fig.  3,  Appendix  1 ). 

5.2  Development  of  near-infrared  center  of  intensity  time  gated  imaging  approach 

A  major  thrust  of  the  project  is  to  develop  optical  imaging  approaches  that  are  simple  and  easy  to 
implement.  During  the  present  reporting  period  we  have  explored  the  development  of  a  new 
near-infrared  optical  imaging  approach  for  locating  a  target  embedded  in  a  turbid  medium.  The 
approach  is  based  on  the  premise  that  a  target  (that  is,  an  optical  inhomogeneity)  within  the 
turbid  medium  alters  the  propagation  of  light  through  the  medium.  Consequently,  the  spatial 
distribution  of  the  output  light  intensity  (SDOLI)  is  different  with  an  embedded  target  than  that 
without  it.  The  target  localization  is  based  on  an  analysis  of  the  spatial  variation  of  the 
transmitted  light  intensity  distribution  for  illumination  at  different  positions  on  the  sample 
boundary.  The  basic  principle,  experimental  arrangement,  and  preliminary  results  are  presented 
in  Appendix  2  (“ Detection  of  a  target  in  a  highly  scattering  turbid  medium  using  near-infrared 
center  of  intensity  time  gated  imaging”),  in  the  form  of  a  manuscript  to  be  submitted  for 
publication. 

The  experimental  arrangement  for  realizing  the  approach  in  practice  is  shown  in  Fig.  2  of 
Appendix  2.  The  scattering  medium  was  a  suspension  of  Intralipid- 10%  in  water,  with  an 
estimated  reduced  scattering  coefficient,  ps  ~  1.168  mm"1  (transport  length  0.86  mm)  and  an 
absorption  coefficient  of  0.0021  mm"1  at  800  nm.  The  Intralipid- 10%  suspension  was  held  in  a 
240  mm  x  160  mm  x  65  mm  rectangular  glass  container.  The  target  was  a  10  mm  x  10  mm  x  3 
mm  neutral  density  filter  (absorption  coefficient  0.23  mm"1  at  800  nm).  Another  sample  was  a 
5  mm  x  5  mm  x  5  mm  piece  of  porcine  liver  (target)  placed  inside  a  150  mm  x  90  mm  x  50  mm 
slab  of  porcine  tissue  held  in  a  rectangular  plastic  cell. 

The  samples  were  illuminated  by  800  nm,  200  ps,  1  kHz  repetition  rate  pulses  from  a 
Ti:sapphire  laser  and  regenerative  amplifier  system.  The  average  beam  power  was  150  mW  and 
the  beam  spot  size  was  approximately  2  mm.  The  laser  beam  was  incident  along  the  z-axis  into 
one  of  the  240  mm  x  160  mm  flat  faces  of  the  container.  The  2 -D  intensity  distribution  of  light 
emergent  from  the  opposite  end  face  of  the  cell  was  recorded  by  an  ultrafast  gated  intensified 
camera  system  (UGICS).  The  UGICS  provides  an  electronic  time  gate  whose  full  width  at  half 
maximum  (FWHM)  duration  can  be  set  as  short  as  80  ps.  The  gate  position  could  be  varied  over 
a  20  ns  range  with  a  minimum  step  size  of  25  ps.  The  sample  cell  was  mounted  on  a  translation 
stage  for  lateral  scanning.  The  input  beam  and  the  UGICS  were  not  scanned. 

The  approach  was  able  to  detect  the  target  in  both  the  cases.  As  shown  in  Fig.  5  of 
Appendix  2,  the  approach  provided  a  two-dimensional  image  of  the  porcine  liver  embedded 
within  the  porcine  tissue. 

6.  KEY  RESEARCH  ACCOMPLISHMENTS 

•  Demonstrated  the  efficacy  of  Optical  Tomography  using  Independent  Component  Analysis 
(OPTICA)  approach  for  detection,  3-D  localization,  and  cross  section  imaging  of  a  tumor 
inside  a  realistic  breast  model  composed  of  excised  breast  tissues  was  detennined  with 
millimeter  accuracy  ( Appendix  1 ). 
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•  Developed  a  new  approach  for  obtaining  two-dimensional  image  of  a  target  embedded  inside 
a  turbid  medium  ( Appendix  2). 

7.  REPORTABLE  OUTCOMES 

Journal  Articles 

1.  M.  Xu,  M.  Alrubaiee,  S.  K.  Gayen  and  R.  R.  Alfano,  “Optical  diffuse  imaging  of  an  ex  vivo 
model  cancerous  human  breast  using  independent  component  analysis,”  IEEE  J.  Select. 
Topics  Quantum  Electron.  14,  43  (2008). 

Presentations 

2.  M.  Alrubaiee,  S.  K.  Gayen,  and  R.  R.  Alfano,  “Near- infrared  time-resolved  and 
spectroscopic  imaging  for  breast  cancer  detection.”  Poster  P42-17  presented  at  the  Era  of 
Hope,  Department  of  Defense  Breast  Cancer  Research  Program  Meeting,  June  25-28,  2008, 
Baltimore,  Maryland.  Abstract  appears  in  p.  293  of  Meeting  Proceedings. 

3.  M.  Xu,  M.  Alrubaiee,  S.  K.  Gayen,  and  R.  R.  Alfano,  “Optical  tomography  using 
independent  component  analysis.”  Poster  P42-17  presented  at  the  Era  of  Hope,  Department 
of  Defense  Breast  Cancer  Research  Program  Meeting,  June  25-28,  2008,  Baltimore, 
Maryland.  Abstract  appears  in  p.  292  of  Meeting  Proceedings. 

8.  CONCLUSION 

The  work  carried  out  during  this  reporting  period  shows  the  potential  for  detection  and  three- 

dimensional  localization  of  a  tumor  within  a  breast  with  significant  accuracy  based  on  the 

differences  in  the  scattering  and  absorption  characteristics  of  the  tumor  and  normal  breast  tissue. 

“ So  What  Section” 

•  A  recent  study  involving  35,3 19  patients  underscores  the  influence  of  primary  tumor  location 
on  breast  cancer  prognosis,  and  makes  it  imperative  that  breast  cancer  detection  modalities 
obtain  three  dimensional  (3-D)  location  of  the  tumor  relative  to  the  axilla.7  The  current  work 
is  an  important  development  in  obtaining  noninvasive  3-D  location  of  a  tumor  within  the 
breast. 

•  Three-dimensional  target  localization  will  enable  closer  probing  of  a  smaller  volume  around 
of  the  target  providing  more  details  since  smaller  pixel  size  could  be  used  without  increasing 
the  computation  time  (as  a  smaller  volume  will  be  probed). 

•  What  is  more  important,  the  work  accomplished  so  far  prepares  us  for  the  next  step,  in  vivo 
optical  breast  imaging  involving  volunteers. 
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Optical  Diffuse  Imaging  of  an  Ex  Vivo  Model 
Cancerous  Human  Breast  Using  Independent 

Component  Analysis 

Min  Xu,  Mohammad  Alrubaiee,  S.  K.  Gayen,  and  R.  R.  Alfano,  Fellow,  IEEE 


Abstract — Optical  imaging  using  independent  component  anal¬ 
ysis  (OPTICA)  has  been  used  for  detection,  3-D  localization,  and 
cross-section  imaging  of  a  tumor  inside  a  model  human  breast 
composed  of  ex  vivo  human  breast  tissues.  OPTICA  uses  a  mul¬ 
tisource  target  illumination  and  multidetector  signal  acquisition 
scheme  to  obtain  multiple  spatial  and  angular  views  of  the  sample 
for  target  localization.  Independent  component  analysis  of  the  per¬ 
turbations  in  the  spatial  light  intensity  distribution  measured  on 
the  sample  boundary  sorts  out  the  signal  originating  from  individ¬ 
ual  targets.  A  back-projection  technique  estimates  the  cross-section 
of  each  target.  The  approach  correctly  provided  the  positions  of  a 
tumor  located  at  the  mid-plane  and  two  glandular  structures  lo¬ 
cated  at  different  positions  within  the  33-mm-thick  model  breast. 
The  reconstructed  cross-section  images  are  in  good  agreement  with 
known  dimensions  of  the  structures,  and  pathological  findings. 

Index  Terms — Breast  cancer,  diffuse  optical  imaging,  indepen¬ 
dent  component  analysis,  near  infrared  (NIR)  imaging,  optical 
mammography,  optical  imaging  using  independent  component 
analysis  (OPTICA). 

I.  Introduction 

NEAR-INFRARED  (NIR)  diffuse  optical  tomography 
(DOT)  is  an  emerging  technology  for  functional  char¬ 
acterization  of  biological  tissues,  and  has  been  actively  investi¬ 
gated  to  image  lesions  in  human  body  organs,  such  as  human 
breast  [1] — [3],  brain  [4] — [7],  and  joints  [8],  [9],  A  state-of-the-art 
DOT  illuminates  the  sample  (consisting  of  targets  embedded  in 
a  turbid  medium)  with  NIR  light,  measures  the  emergent  light  on 
the  boundary  of  the  turbid  medium,  and  uses  an  iterative  image 
reconstruction  method  for  repeatedly  solving  the  forward  model 
of  light  propagation  in  the  medium  with  an  updated  estimation 
of  its  optical  properties  to  match  the  detected  light  intensities. 
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This  problem  of  imaging  targets  in  a  turbid  medium  is  an  ill- 
posed  inverse  problem,  and  a  priori  knowledge  about  the  optical 
properties  of  the  medium  need  to  be  used  to  obtain  a  unique  so¬ 
lution  at  a  cost  of  reduced  resolution  [10]— [13],  Various  prior 
information  such  as  anatomical  structures  obtained  from  X-ray 
or  magnetic  resonance  imaging  (MRI)  and  the  absorption  spec¬ 
tra  of  chromophores  have  been  used  to  improve  the  imaging 
quality  of  the  DOT  [  14] — [16],  The  iterative  image  reconstruc¬ 
tion  is  computation  time  intensive  and  reconstruction  in  2-D 
planar  sections  instead  of  a  3-D  volume  is  commonly  practiced. 
Noniterative  approaches  have  also  been  pursued  [  17] — [19],  Irre¬ 
spective  of  these  developments,  reconstruction  of  images  with 
adequate  spatial  resolution  and  accurate  localization  and  char¬ 
acterization  of  the  targets  remain  a  formidable  task. 

We  have  developed  an  alternative  approach  for  optical  imag¬ 
ing  using  independent  component  analysis  (OPTICA)  [18],  [20] 
that  uses  a  multisource  sample  illumination  and  multidetector 
signal  acquisition  scheme  to  generate  an  extensive  data  set  pro¬ 
viding  a  variety  of  spatial  and  angular  views  of  the  medium.  The 
signals  from  individual  targets  within  the  interrogated  medium 
are  then  sorted  out  by  using  independent  component  analysis 
(ICA)  based  on  their  statistical  independence.  ICA  is  a  statis¬ 
tical  technique  from  information  theory  that  is  able  to  recover 
independent  signals  from  their  measured  mixtures  [21],  [22], 
ICA  has  been  successfully  applied  in  many  biomedical  applica¬ 
tions,  such  as  electroencephalogram  (EEG)  [23]  and  functional 
magnetic  resonance  imaging  (fMRI)  [24],  and  has  been  shown 
to  be  effective  in  separating  signals  from  different  brain  activity 
centers.  In  DOT,  excess  light  absorption  or  scattering  by  the 
individual  targets  embedded  in  the  medium  serve  as  the  source 
of  independent  signals  whose  weighted  mixture  is  recorded  by  a 
detector  on  the  boundary  of  the  medium.  Since  an  independent 
component  originating  from  any  particular  target  relates  directly 
to  how  light  propagates  from  the  source  to  the  target  and  from 
the  target  to  the  detector,  the  recovered  independent  components 
can  serve  as  the  starting  point  for  3-D  localization  and  optical 
characterization  of  individual  targets  in  the  medium.  Such  a 
staged  procedure  has  been  shown  to  significantly  improve  the 
sensitivity  to  small/weak  absorptive,  scattering  and/or  fluores¬ 
cent  targets,  and  can  achieve  a  3-D  localization  of  the  targets 
with  remarkable  accuracy  and  resolution  [18],  [25],  [26]. 

The  independent  component  is  proportional  to  the  strength 
of  the  target  (the  product  of  the  difference  in  the  absorp¬ 
tion/scattering  coefficient  between  the  target  and  the  back¬ 
ground,  and  the  volume  of  the  target)  and  the  convolution  of 
the  light  propagators  from  the  source  to  the  target  and  from  the 
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target  to  the  detector.  The  two  light  propagators  can  be  decon- 
voluted  in  the  Fourier  space.  A  2-D  cross-section  image  of  the 
target  is  obtained  by  back  projecting  the  independent  compo¬ 
nent  onto  the  transversal  plane  at  the  axial  location  of  the  target. 
Every  independent  component  retrieved  by  ICA  represents  the 
signal  from  only  one  target  with  localization  determined  from 
earlier  stage  of  analysis.  So,  a  back  projection  formalism  with 
little  or  no  regularization  can  be  applied  to  obtain  a  cross-section 
image  of  the  target  with  improved  spatial  resolution  than  what 
is  feasible  in  a  conventional  DOT. 

We  have  previously  tested  the  efficacy  of  OPTICA  on  sam¬ 
ples  consisting  of  absorbing  or  scattering  targets  within  tissue 
phantoms  and  fluorescent  targets  in  ex  vivo  tissue  [18],  [25], 
[26].  In  this  paper,  we  use  OPTICA  to  investigate  a  tumor  and 
other  structures  embedded  in  a  “realistic”  model  breast  assem¬ 
bled  using  ex  vivo  human  breast  tissues,  as  a  prelude  to  in  vivo 
breast  imaging.  The  remainder  of  the  paper  is  organized  as  fol¬ 
lows.  Section  II  presents  the  theoretical  formalism  of  OPTICA 
and  the  back-projection  approach  for  obtaining  the  cross-section 
image  of  a  target.  Section  III  describes  the  experimental  arrange¬ 
ment,  method,  and  parameters.  Experimental  results  appear  in 
Section  IV.  The  implications  are  discussed  in  Section  V. 

II.  Theoretical  Formalism  of  Optical  Imaging  Using 
Independent  Component  Analysis 

The  presence  of  targets  (optical  inhomogeneities)  inside  a 
turbid  medium  perturbs  the  spatial  intensity  distribution  of  light 
emergent  from  the  medium  under  illumination  by  a  probing 
beam.  When  illuminated  by  a  point  source  of  unit  power,  the 
change  in  the  light  intensity  distribution  on  the  boundary  of  the 
specimen  due  to  absorptive  and  scattering  targets  can  be  written 
as  [27],  [28] 

-AI(rd,rs)  =  J  d3r  6pa(r)  cG( rd,r)  G( r,rs) 

+ J  d3rSD(r)  cVr  G(rd,r)VrG(r,rs)  (1) 

in  the  first-order  Born  approximation  assuming  that  light  dif¬ 
fuses  inside  the  medium  [29],  Here,  rs  and  rd  are  the  posi¬ 
tions  of  the  source  and  the  detector  on  the  boundary,  5pa  (r)  = 
fia  (r)  —  paa  and  SD(r)  =  D( r)  —  D0  are  the  differences  in  ab¬ 
sorption  coefficient  and  diffusion  coefficient,  respectively,  be¬ 
tween  the  target  at  r  and  the  background  medium,  c  is  the  speed 
of  light  in  the  medium,  and  G(r,r')  is  the  Green’s  function 
describing  light  propagation  from  r'  to  r  inside  the  medium  of 
absorption  coefficient  /./,„(J  and  diffusion  coefficient  D0 . 

OPTICA  assumes  each  inhomogeneity  within  the  turbid 
medium  to  be  a  virtual  source  and  expresses  the  change  of 
the  light  intensity  on  the  boundary  of  the  specimen  as 

- AI(rd,rs )  =  y^  aj(rd)  Sj(rs)  (2) 

3 

where  Sj(rs)  represents  the  jth  target  illuminated  by  the  inci¬ 
dent  wave  at  r s  and  a;  (r,/)  is  the  weighting  matrix  describing  the 
propagation  of  light  from  the  jth  inhomogeneity  to  the  detector 
at  rd.  Each  absorptive  inhomogeneity  contributes  one  term  in 


(2),  and  each  scattering  inhomogeneity  contributes  three  terms 
in  (2)  [18],  The  detected  change  of  the  light  intensity  —A I  is, 
hence,  a  linear  mixture  of  signals  where  at  and  Sj  can  now  be 
interpreted  as  the  jth  weighting  matrix  and  virtual  source,  re¬ 
spectively.  Owing  to  the  statistical  independence  between  these 
virtual  sources,  independent  component  analysis  of  —A I  will 
yield  a  list  of  independent  components  and  recover  both  a,  and 
Sj .  Here,  a,j  and  Sj  are  the  independent  intensity  distribution  on 
the  detector  and  source  planes,  respectively,  for  the  jth  target. 
The  number  of  the  leading  independent  components  gives  the 
number  of  objects.  The  location  of  the  jth  target  is  obtained 
from  the  analysis  of  the  retrieved  independent  component  (Sj 
and  (ij)  that  relates  directly  to  the  source-to-object  and  object- 
to-detector  Green’s  functions  G'fr, .  rsj  and  G(rd, rj  and  the 
optical  property  of  the  target  where  r(  is  the  position  of  the  jth 
object  [18],  [20],  [25],  [26]. 

For  the  slab  geometry  investigated  here,  there  are  three 
virtual  sources  of  specific  patterns  (one  centrosymmetric 
and  two  dumbbell-shaped)  associated  with  each  scattering 
inhomogeneity,  whereas  only  one  centrosymmetric  virtual 
source  is  associated  with  each  absorptive  inhomogeneity. 
Among  the  three  virtual  sources  associated  with  a  scatter¬ 
ing  inhomogeneity,  the  centrosymmetric  virtual  source  is  the 
strongest  and  more  amenable  to  detection  in  a  thick  turbid 
medium  [25].  The  centrosymmetric  virtual  source  and  the 
corresponding  weighting  matrix  are  Sj  oc  G(rj,  rs)  and  aj  oc 
G(rd,rj),  and  Sj  oc  dG/dz(rj,rs)  and  aj  oc  dG/dz{rd,Yj), 
respectively,  for  absorptive  and  scattering  inhomogeneities.  A 
simple  least  square  fitting  of  the  centrosymmetric  component, 
such  as 

minr .  >otj  i/3 .  | [aj1  Sj  {rs )  -  G{ r, ,  rs )] 2 

+  12  [/^iM -G(r<Fri)]2|  (3) 

r  d  J 

for  the  absorptive  object,  can  be  used  to  yield  the  3-D  location 
r  j  and  the  strength  a:l  f3j  of  the  target.  When  a  priori  knowledge 
about  the  property  of  the  target  is  not  available,  (3)  can  still 
be  used  to  estimate  the  3-D  location  of  the  target  regardless  of  the 
absorption  or  scattering  property  of  the  target.  This  is  due  to  the 
fact  that  dG/dz(jCj,rs)  ~  —nG(rj,rs)  and  dG/dz(rd,  r,j  ~ 
—nG(rdlrj),  where  n  =  \/(pa0  —  ico/c)/D0  is  chosen  to  have 
a  nonnegative  real  part  with  ui  the  modulation  frequency  of  the 
incident  wave. 

The  signal  from  the  jth  target  is  simply  given  by  -A/,  = 
aj  (r d)sj  (r , ) .  On  the  other  hand,  the  centrosymmetric  signal  of 
the  jth  target  can  be  approximated  as  a  double  convolution 

-A Ij  (rd ,  rs)  =  J G(pd  -  p,zd,  z-)  Xj(p)  G{p  -  p. ,  zj ,  zs)dp 

(4) 

where  the  integration  is  over  the  z  =  Zj  plane,  Aj  represents  the 
target,  and  pd  and  ps  are  the  lateral  coordinates  of  the  detector 
and  the  source,  respectively.  The  cross-section  image  of  the 
jth  target  Aj  is  a  2-D  distribution  of  the  absorption/scattering 
coefficient  of  the  target  on  the  2  =  Zj  plane.  In  the  Fourier  space, 
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Fig.  1 .  (a)  Schematic  diagram  of  the  experimental  arrangement.  CCD,  charge  coupled  device;  PC,  personal  computer,  (b)  Expanded  view  of  the  sample  input 

(source)  plane  and  exit  (detection)  plane  showing  the  grid  points  in  the  x—y  plane,  (c)  Typical  raw  CCD  image  of  the  detection  plane,  and  how  it  is  cropped  and 
binned  for  analysis. 


the  target  function  Xj  can  be  obtained  from  (4)  as 


At(q)  = 


AI;(q-  qs,q.s) 

G(q  —  qs,  Zd,  Zj ) G*  (qs ,  Zj ,  zs ) 


(5) 


where  q  and  qs  are  the  spatial  frequency  on  the  lateral  plane  and 
denotes  complex  conjugate.  We  choose  qs  =  0  in  the  eval¬ 
uation  of  the  target  function  (5)  since  sources  are  usually  much 
sparser  than  detectors  in  our  setup  where  a  charge-coupled  de¬ 
vice  (CCD)  camera  is  used  to  detect  the  emergent  light  intensity 
on  the  surface  of  the  medium.  The  inverse  Fourier  transforms  of 
Xj  (q)  yields  the  high-resolution  cross-section  image  of  the  /th 
target  due  to  the  high  density  of  detecting  pixels  of  the  CCD.  The 
size  of  the  target  is  estimated  by  the  full- width  at  half-maximum 
(FWHM)  of  the  cross-section  image  Xj . 

To  sum  up,  OPTICA  first  detects  and  retrieves  independent 
components  corresponding  to  each  target  embedded  inside  a 
turbid  medium,  then  obtains  the  3-D  location  and  strength  of  the 
target  from  these  independent  components,  further  reconstructs 
the  cross-section  image  of  the  target  on  the  transversal  plane 
where  the  target  locates,  and  finally,  the  size  and  the  optical 
property  of  the  target  are  estimated. 


III.  Experiment 

The  experimental  arrangement  for  detection  and  localization 
of  the  tumor  in  the  ex  vivo  model  breast  sample  is  shown  in 
Fig.  1(a).  The  model  breast  was  a  70  mm  x  55  mm  x  33  mm 
slab  composed  of  excised  female  human  breast  tissues  provided 
to  us  by  National  Disease  Research  Interchange  under  an  Inter¬ 


nal  Review  Board  approval  at  the  City  College  of  New  York. 
The  model  breast  was  assembled  using  two  pieces  of  ex  vivo 
human  breast  tissues.  The  larger  piece  was  normal  tissue  that 
included  mainly  adipose  tissue  and  streaks  of  fibroglandular  tis¬ 
sues.  The  existence  of  the  fibroglandular  tissues  was  not  known 
prior  to  making  the  measurements. 

The  second  piece  was  mainly  a  tumor  (infiltrating  ductal  car¬ 
cinoma)  with  a  small  amount  of  normal  tissues  in  the  mar¬ 
gins  with  an  overall  approximate  dimension  of  8  mm  x  5  mm  x 
3  mm.  An  incision  was  made  in  the  mid-plane  (along  the  z-axis, 
which  was  the  shorter  dimension  of  the  tissue)  of  the  normal 
piece,  and  some  amount  of  the  normal  tissue  was  removed  from 
the  central  region  making  a  small  pouch.  The  tumor  piece  was 
then  inserted  into  the  pouch,  and  the  incision  was  closed  by 
moderate  compression  of  the  composite  consisting  of  the  nor¬ 
mal  tissue  and  the  tumor  along  .ryz-directions.  The  breast  tissue 
slab  was  contained  inside  a  transparent  plastic  box.  One  of  the 
sides  of  the  box  could  be  moved  to  uniformly  compress  the 
tissue  along  the  z-axis  and  hold  it  in  position.  The  resulting 
specimen,  a  70  mm  x  55  mm  x  33  mm  slab,  was  treated  as  one 
entity  in  the  subsequent  imaging  experiment.  The  position  of 
the  tumor  within  the  slab  was  known  since  it  was  placed  in  the 
position  as  discussed  earlier.  One  of  the  tests  of  the  efficacy  of 
this  imaging  approach  was  to  see  how  well  the  known  position 
is  assessed. 

A  200  yim  optical  fiber  delivered  a  784  nm,  300  mW 
continuous-wave  beam  from  a  diode  laser  for  sample  illumi¬ 
nation.  The  beam  was  collimated  to  a  1  mm  spot  onto  the 
entrance  face  (henceforth  referred  to  as  the  “source  plane”) 
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of  the  slab  sample.  Multiple  source  illumination  was  realized 
in  practice  by  step  scanning  the  slab  sample  across  the  laser 
beam  in  an  xy  array  of  grid  points  using  a  computer-controlled 
translation  stage.  The  xy  array  was  22  x  16  with  a  step  size 
of  2.0  mm.  The  signal  from  the  opposite  face  of  the  sample 
(henceforth  referred  to  as  the  “detection  plane”)  was  collected 
by  a  camera  lens  and  projected  onto  the  sensing  element  of 
a  cooled  16  b,  1024  x  1024  pixel  CCD  camera.  Although  the 
scanned  area  is  42  mm  x  30  mm  on  the  source  plane,  the  im¬ 
aged  area  of  the  detection  plane  was  much  larger,  covering 
the  entire  70  mm  x  55  mm  transverse  area  of  the  model  breast. 
Each  illuminated  pixel  of  the  CCD  camera  could  be  regarded 
as  a  detector.  For  illumination  of  every  scanned  point  on  the 
source  plane,  the  CCD  camera  recorded  the  diffusely  trans¬ 
mitted  2-D  intensity  pattern  on  the  detection  plane.  Each  im¬ 
age  acquisition  took  100  ms,  and  one  stepping  of  the  transla¬ 
tional  stage  took  1  s.  A  total  of  352  images  were  completed 
within  7  min.  The  OPTICA  reconstruction  and  cross-section 
imaging  is  expected  to  be  completed  within  2  min  once  fully 
automated. 

IV.  Results 

A  typical  2-D  raw  image  of  transmitted  light  intensity  distri¬ 
bution  on  the  detector  plane  for  illumination  at  a  typical  scanning 
position  is  shown  in  Fig.  1(c).  The  average  of  all  the  22  x  16  im¬ 
ages  was  used  to  obtain  the  optical  property  of  the  slab  of  breast 
tissue.  The  radial  profile  of  the  intensity  of  the  transmitted  light 
on  the  average  image  was  fitted  to  that  predicted  by  a  diffusion 
model  of  light  propagation  inside  a  slab.  The  transport  mean  free 
path  was  assumed  to  be  1  mm,  the  value  for  a  typical  human 
breast  tissue  at  785  nm.  The  reduced  scattering  coefficient  was 
then  1  mm-1 .  From  the  decay  of  the  radial  profile  of  the  inten¬ 
sity  of  the  transmitted  light,  the  average  absorption  coefficient 
of  the  entire  model  breast  is  found  to  be  //„  =  0.0039  nrnT1. 
Each  raw  image  is  first  cropped  to  retain  the  region  within  the 
window  of  50.4  mm  x  5 1 .3  mm  (out  of  a  total  70  mm  x  55  mm 
transverse  area  of  the  model  breast)  over  which  image  recon¬ 
struction  would  be  performed.  The  size  of  1  pixel  in  the  raw 
image  is  187  /im  x  187  pm.  The  raw  images  are  binned  by 
merging  5x5  pixels  into  one  to  enhance  the  SNR,  resulting 
in  a  total  of  352  images  of  54  x  55  pixels  each.  All  the  binned 
images  corresponding  to  illumination  of  the  grid  points  in  se¬ 
quence  were  then  stacked,  and  used  as  input  for  independent 
component  analysis. 

The  independent  light  intensity  distributions  obtained  by  OP¬ 
TICA  is  displayed  in  Fig.  2(a).  The  3-D  location  of  the  targets 
were  obtained  from  least  squares  fitting  using  (3).  The  fittings 
of  the  independent  light  intensities  over  lines  passing  through 
the  maximum  value  and  along  the  horizontal  direction  are  dis¬ 
played  in  Fig.  2(b).  The  tumor  C  is  found  at  14.8  mm  from 
the  detection  plane  and  centered  at  (33.3,  21.5,  18.2)  mm.  In 
addition,  two  glandular  sites  were  identified.  The  first  glandular 
site  A  is  found  to  be  located  at  2.5  mm  from  the  detection  plane 
and  centered  at  ( 1 1 .2,  22.4,  30.5)  mm;  the  second  glandular  site 
B  is  at  14.6  mm  from  the  detection  plane  and  centered  at  (21.5, 
37.3,  18.4)  mm.  Comparison  of  known  and  3-D  positions  ob¬ 


tained  from  OPTICA  for  the  cancer  site  and  two  glandular  sites 
is  given  in  Table  I. 

The  cross-section  image  of  the  tumor  obtained  from  a 
2-D  inverse  Fourier  transform  of  (5)  is  shown  in  Fig.  3  (left 
pane).  The  right  pane  of  Fig.  3  displays  the  intensity  profiles 
of  the  cross-section  image  along  the  x-  and  //-directions  de¬ 
noted  by  the  white  dashed  lines.  The  FWHM  values  of  the 
intensity  profiles  yield  estimates  of  the  lateral  dimensions  of  the 
tumor  to  be  10.3  mm  x  7.4  mm,  while  the  known  dimensions 
are  8  mm  x  5  mm.  Histological  micrograph  of  the  suspect  site 
confirmed  tumor.  Similar  back-projection  cross-sectional  im¬ 
ages  and  histological  micrographs  were  obtained  (not  shown 
here)  for  the  glandular  tissues  as  well  and  their  transverse  sizes 
were  estimated  from  OPTICA.  The  existence,  location,  and  size 
of  the  glandular  tissues  were  not  known  a  priori.  The  glandular 
structure  A  near  surface  is  estimated  to  be  2.7  and  1.6  mm  in 
size  along  the  x-  and  //-directions  from  the  cross-section  image, 
respectively.  The  size  of  the  glandular  structure  B  at  the  mid¬ 
plane  is  8.7  and  9.2  mm  in  size  along  the  x-  and  //-directions, 
respectively. 

Low  regularization  was  used  in  generating  the  cross-section 
images  in  Fig.  3  to  achieve  maximal  spatial  resolution.  The 
artifacts  in  the  cross-section  images  can  be  suppressed  with  a 
higher  regularization  at  a  cost  of  lower  spatial  resolution.  Since 
the  target  has  been  localized  in  the  earlier  stage  of  analysis,  the 
target  will  not  be  confused  with  artifacts  in  the  cross-section 
images  and  low  regularization  is  beneficial  here. 

The  investigated  ex  vivo  breast  sample  contained  minimal 
amount  of  blood,  and  hence,  the  reconstructed  images  are  for 
the  scattering  property  of  the  sample.  The  change  of  the  re¬ 
duced  scattering  coefficient  p!s  for  the  targets  can  further  be 
estimated  from  the  reconstructed  independent  components  for 
the  sites  A,  B  and  C.  The  value  of  <) p!s  is  given  by  the  ra¬ 
tio  of  the  strength  of  the  target  and  its  volume.  The  sites  A 
and  B  have  lower  scattering  while  the  site  C  has  enhanced 
scattering  compared  to  the  background  (mainly  adipose  tis¬ 
sue).  The  values  of  Sp!s  are  ~0.2  and  0.4  mm-1  for  the 
tumor  and  glandular  tissues,  respectively.  Subsequent  patholog¬ 
ical  analysis  confirmed  the  site  C  as  infiltrating  ductal  carci¬ 
noma,  and  identified  the  other  two  structures  as  glandular  breast 
tissues. 

V.  Discussion 

The  results  of  the  experiments  clearly  demonstrate  that 
OPTICA  can  locate  the  tumor  inside  the  model  breast  with 
high  accuracy.  As  can  be  seen  from  Table  I,  the  lateral  positions 
of  the  tumor  agree  within  0.5  mm,  while  the  axial  position  agree 
within  ~  1  mm  of  the  known  values.  Similar  high  accuracy  in 
the  respective  positions  of  the  two  pieces  of  glandular  tissues  is 
observed  as  well.  The  accuracy  of  the  lateral  positions  does  not 
depend  significantly  on  the  depth  of  the  targets,  while  that  of  the 
axial  position  shows  a  weak  dependence.  For  the  target  located 
close  to  the  detection  plane  (glandular  site  A  at  a  distance  of 
2.5  mm  from  the  detection  plane),  the  axial  position  is  deter¬ 
mined  exactly,  while  for  targets  in  the  midplane  that  are  much 
more  challenging  to  locate,  the  accuracy  is  within  1  mm.  Given 
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Fig.  2.  (a)  Independent  intensity  distribution  on  the  detector  plane  ( z  =  33  mm)  obtained  by  OPTICA  for  the  tumor  C  (left  pane)  and  the  glandular  structures 

A  and  B  (right  pane),  (b)  Corresponding  bottom  panes  show  the  Green’s  function  fits  (solid  lines)  to  the  horizontal  spatial  profile  (denoted  by  circles  and  crosses) 
through  the  center  of  the  intensity  distributions  along  the  dashed  lines. 


TABLE  I 

Comparison  of  Known  and  Optica  Estimated  Target  Locations 


Target 


Known  Position 
(x.  y,  z)  (mm) 


Cancer  Site  (C)  (33,21,16.9) 

Glandular  Site  (A)  (1 1,22,30.5) 

Glandular  Site  (B)  (21,37,17) 


OPTICA 

Estimated  Position 
(x.y.z)  (mm) 

(33.3,21.5,18,2) 

(1 1.2,22.4,30.5) 
(21.5,37.3,18.4) 


that  light  propagation  is  highly  diffusive  in  breast  tissues,  this 
level  of  accuracy  is  quite  significant. 

The  back-projection  formalism  estimates  the  FWHM  values 
of  the  lateral  dimension  of  the  tumor  to  be  10.3  and  7.4  mm 
in  size  along  the  x-  and  //-directions,  respectively,  whereas  the 
known  dimension  is  8  mm  x  5  mm.  This  result  is  expected  due 


to  diffusion  of  light  in  the  tissue,  and  is  in  line  with  the  results 
that  we  obtained  in  our  earlier  OPTICA  studies  [26]. 

Another  important  finding  was  that  OPTICA  predicted  dif¬ 
ferent  scattering  properties  for  the  adipose  tissue  (medium),  the 
tumor,  and  the  glandular  tissues.  The  glandular  tissues  were 
found  to  be  less  scattering  than  the  adipose  tissues  at  the  wave¬ 
length  of  interrogation,  i.e.,  784  nm.  The  tumor  was  found  to 
be  more  scattering.  These  observations  are  consistent  with  the 
known  literature  values  of  scattering  properties  of  different  types 
of  tissues  [30]. 

The  nature  of  the  inhomogeneity  (either  absorptive  or  scat¬ 
tering  or  mixed)  can  be  discerned  by  OPTICA  with  continuous- 
wave  measurement  when  the  SNR  is  high  [20],  [25].  When  the 
SNR  is  not  favorable,  the  recovered  independent  component 
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Fig.  3.  Cross-section  image  of  the  tumor  at  the  z  —  18.2  mm  plane  formed  by  back-projection  (left  pane).  Right  pane:  Spatial  profiles  of  the  cross-section  image 
along  the  x-  and  //-directions  shown  by  the  white  dashed  lines  (right  pane).  The  FWHM  of  the  cancer  site  is  10.3  and  7.4  mm  along  the  x-  and  //-directions, 
respectively. 


will  be  due  to  both  absorption  and  scattering  perturbations  at 
the  site  of  the  inhomogeneity.  The  strength  of  the  target  will 
be  proportional  to  5fj,a  +  k25D  =  Sfia  +  (fj,ao  —  ico/c)5D/D0, 
which  provides  a  way  to  discriminate  between  absorption  and 
scattering  if  measurements  of  multiple  modulation  frequencies 
tu  are  available.  The  capability  of  OPTICA  for  separating  ab¬ 
sorption  from  scattering  inhomogeneities  can  be  significantly 
improved  with  a  time-domain  or  frequency-domain  measure¬ 
ment.  Another  enabling  factor  will  be  carrying  out  multispectral 
OPTICA  studies  for  potential  diagnostic  information. 

OPTICA  can  be  used  for  fluorescent  targets  as  well  [26]. 
The  same  experimental  arrangement  may  be  used,  except  for 
the  introduction  of  filters  to  block  the  excitation  beam  and  to 
transmit  the  fluorescence  light.  What  is  even  more  interesting 
is  that,  a  beam-splitter  and  two  detectors  combination  with  ap¬ 
propriate  filters  may  be  used  to  simultaneously  pursue  absorp¬ 
tion/scattering  OPTICA  and  fluorescence  OPTICA  studies  of 
biological  samples  for  obtaining  coregistered  information  from 
dual  probes. 

OPTICA  is  suited  to  detect  small  objects.  Given  its  ability 
to  identify  low-contrast  small  objects,  the  approach  is  expected 
to  be  especially  useful  for  the  detection  of  breast  and  prostate 
tumors  at  their  early  stages  of  growth. 
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Detection  of  a  target  in  a  highly  scattering  turbid  medium  using  near- 
infrared  center  of  intensity  time  gated  imaging 
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Abstract 

A  near-infrared  optical  imaging  approach  for  locating  a  target  embedded  in  a  turbid  medium 
is  introduced.  The  target  localization  is  based  on  an  analysis  of  the  spatial  variation  of  the 
transmitted  light  intensity  distribution  for  illumination  at  different  positions  on  the  sample 
boundary.  The  approach  is  used  to  detect  and  generate  a  two-dimensional  image  of  an  absorbing 
target  embedded  inside  a  model  scattering  medium  of  thickness  50  times  the  transport  mean  free 
path  of  the  medium,  as  well  as,  a  piece  of  porcine  liver  embedded  in  porcine  tissue. 
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Optical  detection  of  a  target  embedded  in  a  highly-scattering  turbid  medium  is  a  challenging 
problem  with  diverse  potential  applications,  such  as,  imaging  of  a  tumor  in  human  breast,  mines 
in  shallow  coastal  water,  and  airborne  surveillance  through  cloud  or  fog  cover.  Light  multiple 
scattered  by  the  intervening  medium  encroaches  into  the  region  of  geometrical  shadow,  and 
obscures  the  direct  transillumination  image  of  the  target.  A  variety  of  time-resolved,  frequency- 
domain,  and  continuous  wave  methods  for  direct  imaging  and  inverse  reconstruction  of  images 
has  evolved  over  the  years.  '  Most  of  these  approaches  are  difficult  to  implement  because  those 
require  complex  experimental  arrangements,  sophisticated  analytical  schemes  and  time-intensive 
numerical  algorithms  for  generation  of  images. 

In  this  article,  we  introduce  a  transillumination  imaging  approach  for  prompt  target 
detection.  The  approach  is  based  on  the  premise  that  a  target  (that  is,  an  optical  inhomogeneity) 
within  the  turbid  medium  alters  the  propagation  of  light  through  the  medium.  Consequently,  the 
spatial  distribution  of  the  output  light  intensity  (SDOLI)  is  different  with  an  embedded  target 
than  that  without  it. 

The  SDOLI  depends  on  where  and  how  the  turbid  medium  is  illuminated,  as  illustrated 
schematically  in  Fig.  1.  In  the  slab  geometry  of  Fig.  1,  a  narrow  beam  of  light  is  incident  on  the 
source  plane  (z  =  0  plane).  A  fraction  of  the  beam  propagates  diffusely  in  the  z-direction,  and 
emerges  through  the  opposite  side,  i.e.,  the  detection  plane  (z  =  d  plane).  The  SDOLI  at  the 
detection  plane  is  a  two-dimensional  (2 -D)  intensity  distribution  that  can  be  measured  by  a 
charge  coupled  device  (CCD)  camera.  The  black  ellipse  represents  a  target  in  the  medium,  and 
the  gray  spot  at  the  detection  plane  is  the  shadow  of  the  target.  For  a  uniform  scattering  medium 
without  the  target,  the  intensity  distribution  of  the  diffused  light  inside  the  medium  is  symmetric 
about  the  incident  direction,  as  shown  by  the  solid  curve  of  intensity  distribution  in  Fig.  1.  When 


there  is  a  target  with  optical  properties  different  than  that  of  the  turbid  medium,  as  shown  by  the 
black  ellipse,  the  transmitted  light  intensity  distribution  will  be  altered  because  of  scattering 
and/or  absorption  by  the  target.  Consequently,  the  spatial  distribution  of  the  light  intensity  at  the 
detection  plane  will  be  distorted  and  different  with  an  embedded  target  than  that  without  it,  as 
shown  schematically  by  the  dashed  curve.  It  is  possible  to  locate  the  target  within  the  turbid 
medium  through  an  analysis  of  this  intensity  variation.  The  analysis  starts  with  a  calculation  of 
the  center  of  the  SDOLI  at  the  detection  plane  given  by 


M  = 


(1) 


where  /  is  the  light  intensity  recorded  by  a  CCD  camera  and  x  is  the  position.  For  a  uniform 
medium  illuminated  by  a  point  source,  the  transmitted  intensity  distribution  centers  on  the 
incident  direction.  When  the  target  is  located  to  the  right  side  of  the  incident  beam  line,  as  shown 
in  Fig.  1,  it  will  introduce  a  shadow  at  the  right  part  of  the  output  beam  and  lead  to  a  decrease 
and  distortion  in  the  output  intensity  distribution.  This  results  in  the  movement  of  the  intensity 
peak  of  the  output  beam  to  the  left  direction.  Consequently,  the  center  of  intensity  distribution  M 
will  show  a  left  (negative)  movement.  When  the  target  is  to  the  left  of  the  incident  beam,  it  will 
introduce  a  shadow  at  the  left  side  of  the  output  beam  and  M  will  show  a  right  (positive) 
movement.  When  the  object  is  on  the  incident  beam  path,  the  center  of  intensity  of  the  output 
beam  profile  would  be  same  as  that  without  target  inside,  but  the  peak  intensity  would  be 
reduced. 


The  experimental  arrangement  for  realizing  the  approach  in  practice  is  shown  in  Fig.  2.  The 
scattering  medium  was  a  suspension  of  Intralipid- 10%  in  water,  with  an  estimated  reduced 


scattering  coefficient,  jus’~  1.168  mm'1  (transport  length  0.86  mm)  and  an  absorption  coefficient 
of  0.0021  mm  1  at  800  nm.9  The  Intralipid- 10%  suspension  was  held  in  a  240  mm  x  160  mm  x  65 
mm  rectangular  glass  container.  A  10  mm  x  10  mm  x  3  mm  neutral  density  filter  (absorption 
coefficient  0.23  mm'1  at  800  nm)  was  placed  inside  the  medium  as  the  target. 

The  scattering  medium  was  illuminated  by  800  nm,  200  ps,  1  kHz  repetition  rate  pulses 
from  a  Ti:sapphire  laser  and  regenerative  amplifier  system.  The  amplified  pulses  were  not 
compressed.  The  average  beam  power  was  150  mW  and  the  beam  spot  size  was  approximately  2 
mm.  The  laser  beam  was  incident  along  the  z-axis  into  one  of  the  240  mm  x  160  mm  flat  faces  of 
the  container.  The  2 -D  intensity  distribution  of  light  emergent  from  the  opposite  end  face  of  the 
cell  was  recorded  by  an  ultrafast  gated  intensified  camera  system  (UGICS).  The  UGICS  provides 
an  electronic  time  gate  whose  full  width  at  half  maximum  (FWHM)  duration  can  be  set  as  short 
as  80  ps.  The  gate  position  could  be  varied  over  a  20  ns  range  with  a  minimum  step  size  of  25  ps. 
The  sample  cell  was  mounted  on  a  translation  stage  for  lateral  scanning.  The  input  beam  and  the 
UGICS  were  not  scanned.  An  80  mm  linear  region  of  the  sample  was  scanned  in  4  mm  steps 
providing  2 1  data  points. 

Fig.  3  (a)  shows  the  temporal  profile  of  the  tramsmitted  pulse  without  the  target  inside  the 
medium.  The  80  ps  time  gate  was  centered  at  a  time  that  corresponded  to  the  peak  position  of  the 
transmitted  pulse  profile  indicated  by  an  arrow  on  the  profile.  Two-dimensional  images  were 
recorded  at  each  scanning  position,  a  typical  image  being  shown  in  Fig.  3(b).  Fig.  3(c)  shows  the 
horizontal  intensity  profiles  obtained  by  integrating  a  3 -pixel  wide  area  around  the  white  dashed 
line  shown  in  Fig.  3(b).  The  solid  line  in  Fig.  3(c)  is  a  spatial  profile  without  the  target  inside  the 
medium,  while  the  dotted  line  is  the  corresponding  profile  when  the  target  was  on  the  incident 


beam  path.  The  solid  (dashed)  line  in  Fig.  3(d)  presents  the  spatial  profile  when  the  target  was  to 
the  right  (left)  of  the  incident  beam  path. 

For  a  quantitative  evaluation  of  this  intensity  distortion,  the  spatial  intensity  profiles  were 
obtained  for  all  21  scanning  positions,  and  the  corresponding  centers  of  the  intensity 
distributions,  Mt  were  calculated  using  Eq.  (1).  The  average  of  the  21  profiles  was  then  obtained, 
and  its  center  of  intensity  was  calculated  and  taken  to  be  the  reference,  M0.  The  differences,  AM,- 
=  Mi  -  Mo,  were  obtained  and  plotted  as  a  function  of  scanning  position  as  shown  by  the  squares 
in  Fig.  4.  AM  is  negative  in  the  left  side  of  this  center-of-intensity  curve  (CIC),  passes  through  a 
minimum,  reaches  zero,  rises  to  a  positive  maximum  and  tends  towards  zero  further  out  to  the 
right.  The  curve  has  inversion  symmetry  about  the  zero  position.  The  zero  point  shows  the 
position  of  the  embedded  target.  When  the  scanning  beam  is  far  away  from  the  object,  AM  is 
close  to  zero.  The  dashed  curve  in  Fig.  4  shows  the  first  derivative  of  the  center  of  intensity 
Q=d(AM)/dx,  which  peaks  at  the  target  position. 

These  results  show  that  as  the  distance  between  the  input  beam  and  target  increases,  AM 
changes  quickly  at  first.  The  maximum  value  of  Q  occurs  at  the  position  when  the  scanning 
beam  passes  through  the  target.  After  reaching  the  optimal  value,  AM  changes  slowly  towards 
zero.  When  the  target  is  far  away  from  the  input  beam,  its  influence  is  small,  and  the  intensity 
distribution  does  not  change.  The  FWHM  of  Q  vs.  scan  position  curve  in  Fig.  4  is  approximately 
12.6  mm,  which  is  close  to  the  real  width  of  the  target,  10mm. 

The  potential  of  this  approach  for  optical  biomedical  imaging  was  next  explored.  The 
target  was  a  5  mm  x  5  mm  x  5  mm  piece  of  porcine  liver  placed  inside  a  150  mm  x  90  mm  x  50 
mm  slab  of  porcine  tissue  held  in  a  rectangular  plastic  cell.  The  center  of  the  target  was  at  a 
distance  of  25  mm  from  both  the  entrance  and  the  exit  faces  of  the  porcine  tissue.  A  100  mm  x 


60  mm  inner  region  of  the  sample  was  scanned  across  the  incident  laser  beam.  The  sample  cell 
was  scanned  with  a  step  size  of  4  mm  along  the  x  direction  (100  mm  segment).  The  cell  was  then 
stepped  by  3  mm  along  the  y  direction,  and  the  x-scan  was  repeated.  The  scanning  pattern 
continued  until  the  entire  region  was  scanned.  Time-gated  2-D  images  were  recorded  for  each 
scanning  position.  CIC  curves  were  generated  for  every  scan  along  the  x  direction,  and  the  value 
of  Q  was  estimated  for  each  scanning  position.  Plotting  of  Q  vs.  scan  position  ( x ,  y)  generated  a 
2-D  gray-scale  pseudo-color  image  as  shown  in  Fig.  5(a).  The  dark  area  with  high  Q  value 
represents  the  porcine  liver  specimen  that  absorbed  800-nm  light  more  strongly  that  the  porcine 
muscle  tissue.  Fig.  5(b)  shows  a  spatial  profile  of  the  image  was  obtained  by  integrating  a  3-pixel 
wide  area  around  the  white  dashed  line  in  Fig.  5(a)  The  FWHM  of  the  profile  is  10  mm,  which  is 
larger  than  the  target  width  of  5  mm,  but  comparable  to  the  resolution  that  other  optical  imaging 
methods  obtain.5 

In  summary,  a  method  based  on  analysis  of  the  transmitted  light  intensity  distribution  has 
been  introduced  for  obtaining  lateral  position  and  2-D  images  of  an  absorbing  target  in  a  highly 
scattering  turbid  medium  including  biological  tissues. 
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Figure  Captions 

Fig.  1 .  A  schematic  diagram  of  the  diffuse  light  transmission  through  a  turbid  medium  showing 
the  target  and  the  “shadow.”  The  inset  at  the  top  shows  intensity  spatial  intensity 
distribution  without  the  target  (solid  curve)  and  with  the  target  (dashed  curve). 

Fig.  2.  A  schematic  diagram  of  the  experimental  arrangement  (Key:  PS  =  pulse  stretcher,  RGA= 
regenerative  amplifier,  UGICS  =  ultrafast  gated  intensified  camera  system) 

Fig.  3.  (a)  Temporal  profile  of  the  transmitted  pulse.  The  arrow  indicates  the  gate  position. 

(b)  A  typical  time-gated  image  (2-D  intensity  distribution)  of  the  output  plane,  (c)  Spatial 
intensity  profile  without  the  target  (solid  line),  and  with  the  target  (dotted  line),  (d) 
Similar  spatial  profiles  as  in  (c),  except  the  solid  (dotted)  curve  corresponds  to  the  target 
being  located  to  the  right  (left)  of  the  incident  beam  line.  The  arrows  C,  Cl,  and  Cr 
represent  the  peak  intensity  positions  in  the  corresponding  spatial  profiles. 

Fig.  4.  Variation  of  the  center  of  intensity  distribution  ( AM)  and  its  first  derivative  Q  as  a 
function  of  the  scan  position 

Fig.  5.  (a)  Two-dimensional  image  of  a  5  mm  x  5  mm  x  5  mm  piece  of  porcine  liver  embedded 
inside  a  slab  of  porcine  tissue;  (b)  horizontal  spatial  profile  of  the  image. 
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